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Two measures are denned to evaluate the coupling strength of smeared interpolating operators 
to hadronic states at a variety of momenta. Of particular interest is the extent to which strong 
overlap can be obtained with individual high-momentum states. This is vital to exploring hadronic 
structure at high momentum transfers on the lattice and addressing interesting phenomena observed 
experimentally. We consider a novel idea of altering the shape of the smeared operator to match 
the Lorentz contraction of the probability distribution of the high-momentum state, and show a 
reduction in the relative error of the two-point function by employing this technique. Our most 
important finding is that the overlap of the states becomes very sharp in the smearing parameters 
at high momenta and fine tuning is required to ensure strong overlap with these states. 

PACS numbers: TBD 



I. INTRODUCTION 

Lattice QCD has enjoyed great success as a tool for 
first-principles hadron-structure calculations. Early pion 
electromagnetic form factor calculations 0, and nu- 
cleon form factor calculations established the for- 

malism and presented first results establishing the chal- 
lenges ahead for obtaining precision form factors to con- 
front experimental data. Nucleon form factors continue 
to be an active area of research [fjl-fl3| and a compre- 
hensive review of recent form factor calculations can be 
found in flil j and references therein. 

In practice, current lattice calculations were limited 
to a momentum transfer of approximately Q 2 = 3 GeV 2 
due to a challenge of increasing statistical errors. Re- 
cently, calculations of the nucleon and pion form factors 
at Q 2 = 6 GeV 2 have been performed using variational 
techniques [l5j . In this paper we explore very high mo- 
mentum states and propose that, with sufficient opti- 
misation of the smearing parameters alone, momentum 
transfers of the order Q 2 — 10 GeV 2 can be accomplished 
in lattice hadron structure calculations. 

Smearing techniques have seen wide spread use in 
many applications in lattice QCD since first being applied 
to fermion operators [H| . The most notable impacts can 
be found in spectroscopy calculations using variational 
methods p7H23 |. In spite of these successes, there has 
been little in the way of the optimisation of smearing 
parameters for high- momenta states. For low-momenta 
states there is no real need for optimization as the over- 
lap of states is typically slowly varying with the smearing 
parameters. In the following we reveal that this is not 
the case for high-momenta states and finely tuned op- 
timization is very beneficial in accessing these states on 
the lattice. 

Isolation of the ground state at high-momentum is es- 
sential to removing otherwise large and problematic ex- 
cited state contaminations. However, suppression of ex- 
cited states through Euclidean evolution alone encoun- 



ters a rapid onset of statistical noise. We introduce two 
different measures to quantify the coupling of a smeared 
operator to the ground state of a proton relative to the 
near-by excited states, and show how these measures de- 
termine the optimal smeared operator for ground state 
isolation early in Euclidean time. 

We also introduce anisotropy into the smeared op- 
erators in the direction of momentum in an effort to 
improve the coupling to these Lorentz-contracted high- 
momentum states. Our results are complementary to the 
variational techniques of Ref. [l5| in that the optimal set 
of smearings for accessing a variety of momenta can be 
combined to create a correlation matrix providing an ef- 
fective basis for eigenstate isolation. 



II. TWO-POINT FUNCTIONS 

The two-point function of a baryon on the lattice in 
momentum space is given by 



(1) 



where Xi and Xi annihilate and create the baryon respec- 
tively at the sink point x and source point and the 
index i admits various spin-flavor structures for the in- 
terpolators. In the case of the proton, the annihilation 
operator is 



Xi 



(2) 



where u and d represent the spinors for the up and down 
quarks respectively and C is the charge conjugation ma- 
trix. It can be shown that, for positive parity states, 
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where the sum over B represents the ground and excited 
states of the baryon. It is common to average the (1,1) 



and (2,2) elements of the Dirac matrix where the signal 
for positive parity states is large. At zero momentum, 
the Dirac matrix contribution is then 1. The coefficient 
A b provides a measure of the total overlap of Xi &t the 
source and Xi a t the sink with the state B. It is the 
product of the source and sink overlaps which may be 
different if different smearings are used at the source and 
the sink. In this investigation the source will be fixed to 
a point source such that variation in As is proportional 
to the variation in the overlap of Xi which will encounter 
a wide range of different sink smearings. 

Each state decays at a rate proportional to the expo- 
nential of its energy. By evolving forward in Euclidean 
time, excited state contributions die away allowing the 
ground state to be isolated. This is less than ideal for 
the calculation of three-point functions that require ef- 
fective ground state isolation close to the source to avoid 
large Euclidean time evolution and an associated loss of 
signal. It is for this reason that various techniques have 
been implemented for earlier Euclidean-time isolation of 
the ground state. 

When calculating the two-point function, it is possible 
to choose the momentum of the baryon. On the finite 
lattice, momentum is quantised 



P 



2n 



{Px,Py,Pz) 



(4) 



where Nl is the spatial extent of the lattice, a is the 
lattice spacing and p x , p y , p z are integers restricted to 
the range 



(5) 



Due to the construction of the discrete fermion propa- 
gator, momentum input into the two-point function be- 
comes proportional to sin(p), therefore, it is only reason- 
able to consider momentum states where 



< 



Nr 



(6) 



such that the dispersion relation is approximately satis- 
fied. 



III. GAUSSIAN SMEARING 

Gaussian smearing is an iterative procedure applied to 
the source or sink of the two-point function in order to 
improve the relative coupling to the ground state of the 
particle. Consider 



Xi+i(x) = F(x,y) Xl (y). 



with lit 



F (x,y) = (1 - a) &xy 



(7) 



(8) 



where a is a constant, which we set to 0.7. We can in- 
troduce anisotropy to the smearing by introducing a new 
constant a x , which will act only in the x direction, the 
expression for the smearing then becomes, 



F(x, y) = (1 - a a ) 5 xy + ^ (u\{x ~ ax) 6 X - X<y + Ui(x) S x+X ^ + ^ ^ ( u t( x ~ a A) $x-p.,y + u n(. x ) <Wa,v) ( 9 ) 



where a = 0.7 and a and a x are normalised such that 
Aa + 2a x 



G 



IV. MEASURES 



Gusken [16j introduced the measure 



R = 



G 2 {t')e 



+m t 



G 2 (0) 



(10) 



(11) 



for quantifying the ground state isolation of a hadron. By 
taking a point, t , sufficiently late in time such that the 
excited state contributions become negligible, the ground 
state can be evolved back to the source via e +m ° * 



to 



evaluate the fraction of G2(0) it holds. However, with 



sufficient smearing, states can contribute negatively to 
the two-point function, allowing this ratio to exceed 1 
and making it difficult to interpret the results. 

The first measure we introduce follows from this idea 
by determining the deviation of G2 (t) from the ideal two- 
point function of a single ground state. It is similar in 
principle to Gusken's measure, however, it is capable of 
taking into account the presence of states with negative 
coupling to the operator. The measure, Mi is defined 
as, 



Mi 



-1 



t f - U 



-Eo(t-to) 



G 2 (t) 



cm 



(12) 



where G(t) = G(t)/G(t ). The factor -1 makes this 
measure maximal when G(t) is a pure exponential of the 
ground state. The energy Eq is determined from a 4 x 4 



2 



source-sink-smeared variational analysis 23] of the zero 
momentum state with the correct dispersion relation ap- 
plied for finite-momentum states. 

Another common method of extracting coupling effec- 
tiveness is to perform a four parameter, two exponential 
fit on a region close to the source of the two-point func- 
tion, i.e. 

Gat = aie- a2t + 6ie~ M . (13) 

However, this method tends to prove unreliable with the 
parameters varying with the fit window. The method is 
limited by the fact that it can not take into account any 
states with higher energy than the two considered. 

The second measure we introduce works similar to this. 
However, the parameters of the ex pon entials are prede- 
termined by a variational analysis [231 ] . This leads to a 
simple linear fit of known exponentials, i.e. 

Gat = A e- E ° f + X 1 1 + A 2 1 . (14) 

We can then find the proportion of the i-th state in the 
two-point function with the measure 



V. LATTICE DETAILS 

Our calculations are performed on configurations of 
size 32 3 x 64 with a lattice spacing of 0.0907 fm provided 
by the PACS-CS collaboration [24|. These lattices have 
2 + 1 sea quark flavours generated with the Iwasaki gauge 
action [2 a ] and the non-perturbatively improved Clover 
fermion action [26] with the K values for the light quarks 
and the strange quark given by 0.13754 and 0.13640 re- 
spectively, and C'sw = 1-715. This gives a pion mass of 
m n = 389 MeV. 

In order to eliminate any bias caused by smearing in 
the source, we use a single set of propagators generated 
with a point source. All of the smearing is then applied 
to the sink, making the two-point functions smearing de- 
pendent. All momentum will be in the x direction, i.e. 
p y = and p z = in Eq. (JU). 

We use a 4 x 4 correlation matrix to extract our excited 
state masses, constructed from the xi operator with 16, 
35, 100 and 200 sweeps of smearing. We choose to use 
the larger basis in order to ensure that the first three 
eigenstate energies are accurately determined. 

We have verified that no multi-particle states are 
present in the variational analysis by applying the single- 
particle dispersion relation to the zero momentum ef- 
fective state masses to successfully predict the effective 
masses of the same states with non-zero momentum. 

Our error analysis is performed with the second-order 
single-elimination jackknife method. Linear fits are per- 
formed using the normal equations with exact matrix in- 
version where possible and singular value decomposition 
otherwise. 
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FIG. 1. The measure from Eq. JTSJ at p x = in Eq. @. 
Deviation from the ideal two-point function increases by a 
factor of 10 less than 30 sweeps from the ideal smearing level, 
as shown in the inset graph 



VI. RESULTS 
A. Isotropic Smearing 

We first calculate the measure from Eq. (jl"2j) where the 
two-point functions have been normalised 1 time slice af- 
ter the source, with t , = 1 and tt = 6. The two-point 
function is calculated at every sweep of sink smearing be- 
tween 1 and 480, up to an rms radius of 13.68 in lattice 
units. For this particular ensemble, the two-point func- 
tion that shows the highest proportion of ground state 
has 136 sweeps of smearing at the sink, or an rms ra- 
dius of 6.92 lattice units as seen in Fig. [T] Also apparent 
is that the effectiveness of the smearing at isolating the 
ground state is significantly reduced fairly close to the 
optimal amount of smearing. At only 30 sweeps away 
from the ideal number of sweeps, the deviation from the 
ideal two-point function has increased by a factor of 10. 

When we move to p x = 1 in Eq. (U|), which gives mo- 
mentum in the x direction of 427 MeV, the ideal number 
of smearing sweeps reduces by just one sweep to 135 (rms 
radius 6.90 lattice units), as shown in Fig.O This can be 
explained by considering the relativistic 7 factor, which 
is given by the ratio of the relativistic energy momentum 
relation and the ground state mass. The fitted ground 
state mass for the proton is Mp = 1.273(21) GeV, giv- 
ing a relativistic energy of Ep\ p= i = 1.343(23) GeV and 
7 = 1.05. Given that all of the excited states are more 
massive, and therefore exhibit less Lorentz contraction 
than the ground state, it is feasible that there is very 
little difference in the probability distribution between 
this state and the zero momentum state, thus the ideal 
amount of smearing should be very similar to the zero 
momentum state. 

At p x = 3 in Fig. [2J the optimal number of smearing 
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FIG. 2. The measure from Eq. (|12[1 at p x = 1 (left) and p x = 3 (right) in Eq. ((I}. There is little difference between the measure 
at p x — and p x = 1, due to the fact that the probability distributions between the two momentum states are nearly identical. 
At p x = 3, the rms radius of the optimal smearing level is smaller by a factor of 0.85 relative to the p x — state, whereas the 
relativistic 7 factor provides a Lorentz contraction factor of 7 -1 = 0.72 
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FIG. 3. The measure from Eq. (|12[1 at p x — 5 (left) and p x = 7 (right) in Eq. (j4|. The value of the measure at the optimum 
number of smearing sweeps for this momentum state is approximately equal to that of the p x = 3 state, indicating that good 
ground state isolation is possible even at higher momenta. At p x = 7, the deviation from the ideal two-point function has 
increased by a factor of 10 only 5 sweeps from the optimal smearing level, as shown in the inset graph. 



sweeps has decreased to 98. The maximum value of the 
measure has also decreased relative to the lower momen- 
tum states, indicating relatively more excited state con- 
tamination, though still achieving good isolation. The 
ratio of the rms radius of the optimal smearing for this 
state to the optimal smearing for the ground state is 0.85, 
compared to the relativistic r y~ 1 factor of 0.72. At p x = 5, 
corresponding to a momentum transfer of approximately 
4.55 GeV 2 , shown in Fig.|3l the optimal number of sweeps 
is 52 (rms radius 4.27 lattice units). However, the maxi- 
mum value of the measure is close to the maximum value 
for the p x — 3 case, indicating that very efficient isolation 
is possible, even at larger momentum transfers. 

Moving to p x = 7, equivalent to a momentum trans- 
fer of 8.93 GeV 2 , there is significant noise far from the 



source in the two-point function, even for highly opti- 
mised smearing values. Hence we consider tt — 5 in the 
measure from Eq. (1121) at this value of momentum. The 
ideal number of sweeps decreases to 27 sweeps, or 3.08 
lattice units rms radius, seen in Fig. [3l Notably, the de- 
viation from the ideal two-point function increases by a 
factor of 10 only 5 sweeps from this optimal value, corre- 
sponding to a change in rms radius of less than 0.3 lattice 
units. 

Using the measure described in Eq. (fl"5|) . we first 
consider the three exponential fit between time slices 
1 and 6 after the source with masses 1.273(21) GeV, 
2.301(28) GeV and 2.786(95) GeV as determined in our 
correlation matrix analysis. From the results in Fig. 21 
we can see that, in the region where the first measure 
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FIG. 4. Ground state proportion from the three exponential 
fit at p x — in Eq. (j4}. There is insufficient information on 
the second excited state close to the optimal amount of smear- 
ing, thus requiring use of the two exponential fit to determine 
the optimal amount of smearing with this measure. 



FIG. 5. Ground state proportion at p x — in Eq. (|4]). Con- 
tamination due to excited states increases rapidly away from 
the optimal smearing level. There is good agreement between 
the two exponential fit here and the three exponential fit in 
Fig. [4] away from the optimum smearing levels. 



predicts ideal smearing levels, there is a sharp change in 
the structure of the graph. In order to determine the 
cause of this, we compare with the fits containing only 
the ground and first excited states. Fig. [5] shows that 
the optimal number of smearing sweeps lies close to the 
value predicted by the first measure. The overlap at the 
optimal number of sweeps, 138 in this case, is 99.31(8)%, 
indicating that, in the three exponential fit, we are at- 
tempting to fit two quickly decaying exponentials using 
only 0.69% of the signal available. This leads us to be- 
lieve that, in the regions of ground state dominance where 
we are most interested, the coefficient from the quickly 
decaying third state cannot be determined accurately, 
therefore dominates well beyond where it should be al- 
lowed to contribute at all. For this reason, we will only 
consider fits using the ground and first excited states. 

The contamination due to excited states in the two ex- 
ponential fit at zero momentum increases rapidly away 
from the optimum smearing level. Of the smearing 
sweeps used to extract the masses from the variational 
analysis, the one that shows the most overlap with the 
ground state is 200 sweeps, or an rms radius of 8.55 lat- 
tice units, with 77.69(7)%, or 32 times more excited state 
contamination than the optimal smearing level. 

At the first non-zero momentum state, the results 
present similarly to the first measure, the optimal amount 
of smearing is 1 sweep less than that of the non-zero mo- 
mentum ground state, and 2 sweeps more than the opti- 
mal amount determined by the first measure. At p x = 3 
in Eq. (U) shown in Fig. [51 the overlap is maximised at 101 
sweeps of smearing, or an rms radius of 5.95 lattice units, 
once again agreeing within only a few sweeps of the op- 
timum level suggested by the first measure. Remarkably, 
considering the use of a point source, the proportion of 
ground state present at this optimal amount of smearing 
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FIG. 6. Ground state proportion at p x — 3 in Eq. Q. As 
momentum increases, the contamination due to excited states 
increases more rapidly away from the ideal smearing level. 



is 98.87(12)%. 

At p x = 5 and p x = 7 in Fig. [7] there is again good 
agreement between the two measures, with the optimal 
smearing level being 53 and 26 sweeps respectively. Even 
at a momentum transfer of 8.93 GeV 2 , 97.20(20)% over- 
lap is achieved with the ground state, and once again, 
very few sweeps from the optimum level, the overlap 
drops dramatically. At p x = 7, far from the optimal 
number of smearing sweeps, it is unlikely that any highly 
Lorentz contracted state would couple to such a large 
sink. The second peak in Fig. [7] can therefore be consid- 
ered to signify a limit to the domain of validity of the 
measure. 
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FIG. 7. Ground state proportion at p x = 5 (left) and p x = 7 (right) in Eq. (j4]) Even at these very high momentum transfers, 
good overlap with the ground state is achieved for an optimised sink. Far from the optimal number of smearing sweeps at 
p x = 7, it is clear that the measure is no longer applicable, as there would be little, if any highly Lorentz contracted ground 
state present. 



B. Anisotropic Smearing 

As anisotropy is introduced to the smearing as de- 
scribed in Eq. ([9]), we consider the first measure from 
Eq. (|12p at the first non-zero momentum state and find 
that there is no improvement to the ground state iso- 
lation, as shown in Fig. [5] There is, however, an ideal 
number of sweeps that increases for decreasing a x that 
shows approximately equal ground-state proportion rel- 
ative to the isotropic smearing case. 

At p x = 3 in Eq. in spite of the clear difference in 
the smearing sweeps required to maximise overlap with 
the source, Fig. |9] shows that introducing anisotropy to 
the smearing does not result in improved isolation of the 
ground state. The structure of the curve is similar to 
that of the p x = 1 state, where there is an optimal num- 
ber of sweeps for every value of a x which increases with 
decreasing a x . 

Once again, there is no improvement in the ability of 
anisotropic smearing to isolate the ground state at the 
momentum of p x = 5, as shown in Fig. 1101 The structure 
revealed in the lower momentum states persists for this 
state and for the p x = 7 state in Fig. QTJ From these 
results, optimisation of the number of smearing sweeps 
alone is sufficient to achieve good isolation of the ground 
state of the two-point function at a range of momenta. 

We now investigate how anisotropic smearing affects 
the signal-to-noise ratio or quality of the two-point func- 
tion at high momenta. Since we have ensured that the 
ground state is isolated as close to the source as possible, 
we now determine the quality of the signal a few time 
slices away from the source. We consider the relative er- 
ror of the two-point function four times slices after the 
source at the optimal number of smearing sweeps for each 
value of our anisotropy parameter, a x . 

For p x = 3, Fig. IT21 shows the two-point function at t = 



4. The smallest relative error occurs when the smearing 
is isotropic. Increasing the momentum to p x = 5 lattice 
units shows that there is only a small improvement to 
the relative error for values of a x ~ 0.48. It is worth 
noting that the first of the minima visible in Fig. (TS] at 
a x — 0.36 corresponds to the anisotropy expected due to 
Lorentz contraction as a x /a = 0.51 equals 7 _1 = 0.51. 

The banding structure visible in Fig. [T3] is a result of 
the optimal number of smearing sweeps increasing for de- 
creasing values of a x . Each discontinuity in the graph for 
a x > 0.36 is the result of the optimal number of smearing 
sweeps decreasing by 1. It is an artifact resulting from 
the density of the points in a x being much finer than the 
density of the points in the number of smearing sweeps. 

Moving to p x — 7 in Fig. [13] we see a distinct im- 
provement in the correlation-function relative error when 
anisotropy is introduced. Both a x = 0.26 and 0.32 pro- 
vide a 10% reduction in the error relative to that observed 
at the isotropic value of 0.7. The values of a x ~ 0.26 to 
0.32 provide a x /a = 0.37 to 0.46, in accord with the 
value of 7 _1 = 0.39 predicted by Lorentz contraction. 



VII. CONCLUSION 

We have presented two new measures of the effective- 
ness of smeared operators in isolating the ground state of 
a hadron in the two-point function. Both measures show 
good agreement with each other. We have performed 
a detailed analysis of ground state isolation with each 
measure and have shown that optimisation of the smear- 
ing can lead to remarkable improvement to the ground 
state isolation. Furthermore, the ability to isolate the 
ground state decreases dramatically a few sweeps from 
the optimal number of smearing sweeps for the higher 
momentum states. In selecting a basis for a correlation 
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matrix analysis, these optimal smearing parameters are 
preferred. 

On the introduction of anisotropy to the smearing, 
we found that there was no appreciable improvement 
to the overlap with the ground state. The relative pro- 
portion of the ground state for an isotropic source is al- 
ready high. Optimising the number of sweeps of isotropic 
smearing alone is sufficient to ensure maximal isolation 
of high-momentum ground states. The introduction of 
anisotropy does provide a small improvement to the cor- 
relation function of high-momentum states a few Eu- 
clidean time slices after the source. 

Our results indicate that future studies of high- 



momentum states should adopt this relatively cheap pro- 
gram of tuning the smearing parameters to optimize iso- 
lation and overlap with the states of interest. We antici- 
pate this approach will be of significant benefit in future 
form factor studies. 
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FIG. 8. The first measure from Eq. (1121) (left) and the Ground State Proportion (right) with anisotropic smearing at p x = 1 
from Eq. (j4}. Introducing anisotropy to the smearing does not improve the isolation of this state. However, the Lorentz 
contraction is small so little improvement would be expected. 




FIG. 9. The first measure from Eq. (|12|) (left) and the Ground State Proportion (right) with anisotropic smearing at p x = 3 
from Eq. Q. No improvement is seen in the isolation of the ground state, in spite of the relativistic 7 factor of 1.39 giving a 
length contraction factor of 0.72 in the x direction. 




FIG. 10. The first measure from Eq. (|12f) (left) and the Ground State Proportion (right) with anisotropic smearing at p x = 5 
from Eq. l|4|. The structure observed in the plots of the p x — 3 state is retained, with more sweeps of smearing required as 
anisotropy is increased. 
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FIG. 11. The first measure from Eq. {12} (left) and the Ground State Proportion (right) with anisotropic smearing at p x = 7 
from Eq. (|4j|. Even at a momentum of 2.99 GeV, anisotropy in the smearing does not improve isolation of the ground state. 
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FIG. 12. Relative error in the two-point function measured 
four time slices after the source for p x = 3 as in Eq. Q. At 
this momentum, isotropic smearing provides the best relative 
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FIG. 13. Relative error in the two-point function measured four time slices after the source for p x = 5 (left) and p x — 7 (right) 
as in Eq. Q. At p x = 5, there is a small amount of improvement for anisotropic smearing at a x /a in the region of ■y~ 1 = 0.51. 
At p x = 7, a 10% improvement in the relative error is seen for values of a x ~ 0.26 to 0.32 where a x /a = 0.37 to 0.46, in accord 
with the value of 7 -1 = 0.39 predicted by Lorentz contraction. Note that the emergent banding structure reflects a change in 
the optimal number of smearing sweeps by one. 
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